A universal set of gates for defect qubits in diamond and silicon carbide 
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We propose a fast scalable optical design for single and two-qubit gates for defect qubits in diamond 
(NV centers) and in silicon carbide, which are promising candidates for room temperature quantum 
computing. Two-qubit gates are mediated by microcavity photons that are coupled to optical 
transitions in these systems. Single-qubit operations are based on optically-assisted selectivity of 
microwave control. This approach does not require tunability of electronic transitions of the defects, 
and the cavity mode frequency remains off-resonance with the optical transitions at all times. The 
gates are controlled by near-resonant narrow-bandwidth optical pulses. We perform full quantum 
numerical modeling of the proposed gates and show that high fidelities can be obtained. 



Qubits encoded by electron states of defects in dia- 
mond and in silicon carbide have become promising can- 
didates for room temperature quantum information pro- 
cessing. Long coherence times of ^ 50 — 200/xs, initial- 
ization, readout and single qubit operations have been 
demonstrated [1] [5] for qubits in both systems at room 
temperature. Simple entangling operations have been 
demonstrated on the system of a negatively charged ni- 
trogen vacancy center (NV) in diamond coupled to a 
nearest carbon (C) [3l |4] or a nearest nitrogen (N) [5j [6] 
nuclear spin. However these systems are not scalable. 
Entangling operations of two distant NV qubits have re- 
cently been performed in a challenging, state-of-the-art 
experiment by joint measurement of the photons emit- 
ted by the two NV centers [7\ . The challenges associated 
with using this type of entanglement in a quantum in- 
formation processor pose stringent requirements on the 
frequencies of the NV centers. Further, the approach is of 
probabilistic nature, — only ~ 3% of the emitted photons 
come from the zero phonon line of the NV center ^ — 
resulting in a low rate of successful operations. It has also 
been proposed that nitrogen defect nuclear spins in dia- 
mond might potentially be used to mediate a long-range 
interaction between distant NV qubits [S]. However, in- 
troduction of a large number of defect nuclear spins as 
well as problems caused by imprecise defect positioning 
in such systems make them impractical. Thus, exper- 
imentally viable deterministic two-qubit gates between 
distant defect qubits remains an important challenge. 

Recent developments in photonic microcavities in both 
diamond [51121 and in silicon carbide lOJ have opened op- 
portunities to couple distant defect qubits via cavity pho- 
tons. Microcavities in these systems have been fabricated 
in either ring 8J or void-pillar geometries [9| with defect 
qubits placed near the surface. High quality factors of 
optical modes and large values of photon coupling to the 
defect optical transitions have been reported for both di- 
amond and silicon carbide systems [H [TOl Hi] . This has 
made photon-mediated distant qubit-qubit interactions 
experimentally accessible. 

We propose a robust optical design for a universal set 
of deterministic gates in diamond and in silicon carbide 



systems using microcavity photons to mediate interac- 
tions. We show that two distant defect qubits can be 
coupled to cavity photons via a two stage optical activa- 
tion scheme. This approach does not require fine tuning 
between cavity mode frequency and optical transition fre- 
quencies of the defects, which remain off-resonance with 
the cavity mode. This leads to effective isolation of the 
single-defect excitations in each qubit, which are used 
to perform fast optical single-qubit rotations and read- 
out. The excited states involving two defects decouple 
from the cavity at substantially larger detunings. These 
states can be used to mediate a qubit-qubit interaction 
and perform a two-qubit control-Z (CZ) gate. By com- 
bining these optical operations with simple global mi- 
crowave rotations we construct a universal set of gates for 
scalable architectures aimed at room temperature quan- 
tum computing. We perform full quantum mechanical 
simulations and find that high values of fidelity can be 
achieved for experimentally reasonable parameters |1H - 

For deflniteness we will focus mainly on the case of 
NV in diamond. The properties of defect centers in 
silicon carbide are similar [Ml [TBHIS] . The NV defect 
has eight states of interest [TH jTS] (see Fig. [l]), six of 
which, denoted from symmetry considerations as 'i^2m. 
and ^Ejn (m — 0,±1), participate in optical transitions. 
The other two ^Ai and typically involve non-radiative 
recombination with emission of phonons 14J. An ex- 
ternal magnetic field, B, splits the states with different 
(to = ±1) spin. The qubit is typically encoded by states 
|0) = |'M2m=o) ^'^d |1) = \^A2m=-i) and is manipulated 
with microwave pulses dlllElE]- The excited states 3£'„ 
are typically used to initialize the qubit to the ground 
state |0) and to read its state. These procedures rely on 
the differences in recombination times for states with dif- 
ferent spin TO and are performed similarly for diamond 
and silicon carbide defect qubits [H [51 [TH [T^ . 

The optical transitions 'M2,„ ^ ^E„i are spin con- 
serving. Their frequencies depend on the local envi- 
ronment and therefore are different for different defects 
[l8l [201 |2T] . Often, in order to couple two qubits via 
a cavity mode, transitions in each qubit have to be 
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Figure 1: Nitrogen- vacancy center (NV) in diamond as a 
qubit. (a) The relevant states of the NV in diamond [14] . 
The qubit is encoded by states ?42m=o ^.uA M2m=-i- Non- 
radiative decay of excited states takes place via and Mi 
(red) [14j. (b) Two such qubit systems with unequal energies 
of excited states '^E interacting via a microcavity photon of 
frequency u)c are used to perform control-Z (CZ) gate. 



tuned dynamically into resonance with the cavity mode 
to enable the interaction for a period of time thus per- 
forming the gate [22l |23]. Instead, we use the natu- 
ral spectral differences between the defects to our ad- 
vantage to perform single- and two-qubit phase-control 
gates via optical pulses. The qubits remain ofF-resonance 
with the cavity. A single qubit phase gate is defined as 
Z{0) = cxp{idaz/2). It acts to add a phase 9 to one of 
the qubit states A\0) + B\l) A\0) + ^11)6*". The sim- 
plest two-qubit phase-control gate, the control-Z (CZ) 
gate, is defined by ^|00) + B\01) + C\10) + D\ll) 
A\00) + B\01) + C\10) - D\n). 

Defects in diamond or silicon carbide in a microcav- 
ity are about a fim away from each other or closer. 
Such proximity makes it difficult to perform individual 
single-qubit operations with microwave pulses. However, 
we show that optically-assisted selectivity of microwave 
control can be achieved even when the microwave field 
is applied to the entire system, performing simultane- 
ous operations on all qubits. The selectivity can be 
achieved by inserting a single-qubit optical Zi{9) gate 
acting on ith qubit between the two conjugate microwave 
operations X and X = X^^ acting on all qubits si- 
multaneously. The particular choice of X is not im- 
portant, and we will use X — X{Tr/2) (S) X{'k/2) and 
X = X{~tt/2) ® X{-tt/2) with xle) = exp(i6'cr^/2) to 
illustrate the procedure. When both conjugate opera- 
tions are performed on a two-qubit state, it remains un- 
affected, i.e., XX I*) = 1*). However, if an optical Zi{9) 
is performed on one (i-th) qubit between X and X, an ar- 
bitrary single-qubit rotation is carried out. For example, 
Y,{9) = y.Z,{9)X and X,{9) = Y,{~n /2)Zi{9)Y,{n /2) 
together with Zi{9) provide an arbitrary single qubit 
gate for the «-th qubit. Single qubit operations together 
with the CNOT = Y2{-tt/2)GZ Y^ij^ji) entangling gate, 
where Y(9) = exTp{i9ay/2), form a standard universal set 
of gates [2^. 

To see how Zi{9) and CZ can be constructed in our 
physical system, we consider the spectrum of the two NV 
defects in a microcavity (see Fig. [2| . In the rotating wave 
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Figure 2: The spectrum of two NV centers interacting with 
a cavity mode, (a) Schematic depiction of the spectrum for 
varying values of uuc- The spectrum can be classified ac- 
cording to the total number of optical excitations. The zero- 
excitations part, shown in panel (b), is used to encode qubits. 
Panel (c) shows the part of the spectrum with one excitation, 
and panel (d) the part with two excitations. Standard sym- 
metry notations are used to label the energy states. In ?42&ci 
each defect is in state and there is one photon in a cavity 
(ci). Notation ^E{1) is used to show that the first defect is 
in state ^E (while the second is in state M2 with no optical 
excitation, and the cavity has no photons). The black dashed 
curves in (c) and (d) mark the excited states used to perform 
CZ gate. 



approximation, the spectrum is given by the Hamiltonian 

<S,i 

^J2{\'M^)m){'E{tUa^ + h.c.) , (I) 



a,i,m=0,±l 



-I- wco + 7 

1=1,2; m=0,±l 

where 2t € {^A2,^E}, <B G {^Ai,^E}, e'^^{i) = \m\AA + 
mgB and e^(z) = A + (5,,2'5A + \m\AE + mgB. The 
values for the energies e<s{i) of the states ^Ai and 
are not important for the following discussion. The 
Hamiltonian ([T]) conserves the total number of excita- 
tions —J2im \^E{i)m){^E{i)m\ + a'^ a and can be diag- 
onalized in each excitation subspace independently. The 
general structure of the states as functions of the cav- 
ity mode frequency ujc is shown in Fig. [2]ja) . The qubit 
subspace states |00), |01), |10) and |11) [Fig.[2];b)] are un- 
affected by the cavity mode. At higher energies electron 
and cavity states mix. The resulting structure shown in 
Fig. [2jc-d) is formed by a collection of simple avoided- 
crossings that involve only a few states due to optical 
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selection rules of the Hamiltonian ([T]). 

In the one-excitation part of the spectrum [Fig. ^c)] 
the states are coupled in groups of three involving two 
anti-crossings: the excited states of each of the two NV's 
in the same spin configuration interact with the cavity 
photon state. When the energies of the excited states 
are sufficiently far from each other, the two anti-crossings 
effectively separate resulting in local (i.e., isolated from 
one another) electron-photon states, e.g., |e) in Fig.[2]Jc). 
In this case the interaction between the two defects de- 
creases rapidly as ~ jAa/SA'^ and for typical values of 
7/(5A ~ Aa/SA ^ 0.01 can be neglected. The decrease 
of the interaction is seen by noting that there are two 
off-resonant transitions involving a cavity photon, each 
adding a factor of 1/ SA. As a result, the single-excitation 
sector is ideal to perform local (single-qubit) optical Z{9) 
operations, as in Ref. ^24j for Zi{0) gate. The transition 
between states |00) and |eO) [Fig. [2]^c)] is addressed with 
a single off-resonant 2tt (or cyclic) pulse. For the above 
values of j/SA the energy difference between states |01) 
and |el) is nearly identical to that in between |00) and 
|eO), and thus only the first qubit is affected by the con- 
trol pulse. An optical pulse field can be added to the 
Hamiltonian as 

np{t-tp)cosujpt J2 {\%i^m)CE{i)m\+h.c.) , (2) 

m=0,±l 

where ftp and ujp are the Rabi and carrier frequencies 
of the laser pulse. For certain shapes of flp{t) [5S] a 
cyclic pulse is produced in which the entire population 
of the state |0) of the first qubit is transfered to state 
|e) and back. The phase accumulated during the pulse, 
|0) — > |e) — |0)e*^, depends on the detuning of ojp from 
the frequency of the transition and on the pulse shape 
rip. A Z2{0) operation is similar and is based on state 
|0e); see Fig.^c). 

We now consider the cavity-induced effective interac- 
tion between the defect qubits in the two-excitation sub- 
space [see Fig. [2|d)]. The effective interaction in the 
two-excitation subspace is greater than that in the one- 
excitation subspace. The energy of the state |"ee"), 
which involves excitations in both defects and a pho- 
ton [see Fig.[2j[d)], differs from the combined energies of 
states |0e) and |eO), which involve excitation of only one 
of the defects. This is the basis for "on-demand" qubit- 
qubit coupling. The state | "ee" ) , however, is not directly 
accessible with an optical field. It can be reached only if 
the system is first transfered into the one-excitation sub- 
space, e.g. |eO). This can be done with a J dtil{t) = tt 
resonant pulse tuned to the corresponding transition. Af- 
ter that, a / dtil{t) — 2tt resonant cyclic pulse tuned to 
the transition between |eO) and | "ee") can be used to per- 
form a Z2{2tt) operation if and only if the first qubit was 
in state |0). Such conditionality is rooted in the difference 
between the energy of | "ee" ) and the sum of energies of 
|0e) and |eO) for non-zero 7. If 7 = this energy differ- 



ence is zero and the state |"ee") is the product |e)|e) and 
the second pulse produces a (unconditional) single-qubit 
Z2{2tt) rotation. After the conditional Z2{2tt) rotation, 
the system must be restored to the qubit subspace with 
the pulse identical to the first one. The pulse sequence 
to perform CZ gate is shown in Fig. [3](a). 

The coherence of the proposed single- and two-qubit 
gates is determined by the fidelity of their components — 
microwave rotations and optical pulse-controlled phase 
gates. The fidelity of the microwave rotations is lim- 
ited primarily by decoherence of the qubits. In diamond 
and silicon carbide systems the time scale of the qubit 
decoherence is large ~ /Lts [H |2]- In the optical gates, 
decoherence can take place at a much shorter time scale. 
In this case there are two sources of possible fidelity loss: 
(i) destructive interference due to transitions that were 
not accounted for, and (ii) real losses due to radiative 
and non-radiative recombination of the excited electron 
states as well as loss of photons from the cavity. Deco- 
herence during a single qubit Zi{9) operation is similar 
to that due to the first pulse in the CZ gate. Therefore a 
lower bound to fidelity of the proposed gates can be ob- 
tained by calculating the fidelity of a single CZ gate. We 
include all eight states shown in Fig. [ija) for each qubit 
and ten states to represent the cavity mode in our numer- 
ical modeling. The fidelity of the CZ gate is evaluated 
from the reduced density matrix of the two-qubit system 
before, p(0), and after, p = p{tg), the gate. The evolution 
is computed using the Bloch-Redfield Master equation 
tp=[H + V{t), p] + J2s i^sLs{p} where = [PspPj ~ 
{PlP,p + pPlPs)/2] dlUS]. We take the state p(0) to 
be a pure state |V'o)(V'o|- The fidelity of the gate opera- 
tion is given by F('(/'o, p{V'o}) = {{ipopcz pUcz\^o)Y''^ 
In quantum computation the state V'o can vary depend- 
ing on the input data and the algorithm involved. To 
estimate the decoherence introduced by the gate in an 
arbitrary circuit, we compute the fidelity averaged over 
the initial state [5U] 

F'^Y. {n\U},,pm{j\}Ucz\m) (3) 

ijnm— {1,4} 

The resulting fidelity is shown in Fig.[3]^b-c). In Fig.[3][b) 
we plot the fidelity optimized over the values of pulse 
bandwidth cr as a function of pulse separation time 5t for 
cavity quality factors in the range A/"fQ = {0.003,0.1} 
p,eV . The values for recombination rates in NV defects 
(~ 10ns) were taken from [M]. We used 7 15 /^eV, 
5 A = 100 ^eV, and gB ^ Q.l peV in Fig. [3| We find 
that the maximum fidelity is reached when pulses over- 
lap thus reducing decoherence. The fidelity as a function 
of a and Q is shown in Fig. [3f[c) for aSt = 3. The fi- 
delity peaks at cr/7 ~ 0.2 and is smaller for higher and 
lower values of a due to unintended dynamics and real 
losses, respectively. The fidelity of the CZ gate applied 
to silicon-carbide defect qubits is qualitatively similar but 
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Figure 3: Fidelity of the CZ gate, (a) Pulse sequence used, 
(b) Variation of fidelity as function of spacing St = t2 ~ ti = 
ta — t2 between the pulses. The maximum value of fidelity 
over the range of ct's and Q's of panel (c) is shown, (c) Fi- 
delity as a function pulse bandwidth, a, and cavity quality 
factor, Q. Each point corresponds to the optimal placement 
of the cavity mode frequency coc relative to the energies of 
the excited states, as shown in panel (d). Inset in panel (d) 
shows fidelity as a function of ujc at A/jQ = 1.3 x 10"'^ and 
with (t/7 — 0.17,0.05, and 0.03, given by solid, dashed, and 
dashed-doted curves respectively. 



involves recombination rates specific to that system. In 
both figures we chose optimal values of the cavity fre- 
quency uJc with respect to the excitation energies of the 
defects. The optimal ujc for each point of Fig. [sjc) is 
shown in Fig. [sj^d). Typically the gate fidelity as a func- 
tion of uc has a wide plateau of high values [see the inset 
in Fig. [sj^d)], and then it decreases when ujc is far above 
or below A. As a result, precise positioning of the defect 
and cavity energy levels with respect to each other is not 
necessary. A similar procedure for a CZ gate can be for- 
mulated using the second qubit as a conditional qubit, in 
which case high fidelities will be reached at luc — A < 
in a symmetric fashion. In both cases the fidelity of the 
gate becomes small when ujc is far off resonance with A 
and 7/(5 A becomes vanishingly small. The maximum fi- 
delity of the gates is limited primarily by the cavity Q. 
Thus, in order to operate, the gates require the strong 
coupling regime, which is given by A/7Q ^ 1. 

We have demonstrated that the fidelity of the proposed 
universal set of gates can reach high values. The single 
and two-qubit operations are performed using different 
classes of states in the excitation spectrum of the system, 
and are not very sensitive to the positioning of the cav- 
ity mode frequency with respect to the defect energies. 
This can be used to couple multiple spectrally distinct 
defects pairwise using a single cavity mode without loss 
of fidelity due to their cross-talk. Thus, this cavity-based 
design has an advantage of scalability in these systems. 



Additional scalability can be obtained by coupling sev- 
eral cavity systems via wave guides. 
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